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ABSTRACT 

The RXTE satellite observed the Coma cluster for ~177ks during 
November and December 2000, a second observation motivated by the 
intriguing results from the first ~87ks observation in 1996. Analysis 
of the new dataset confirms that thermal emission from isothermal gas 
does not provide a good fit to the spectral distribution of the emission 
from the inner 1" radial region. While the observed spectrum may be 
fit by emission from gas with a substantial temperature gradient, it 
is more likely that the emission includes also a secondary non-thermal 
component. If so, non-thermal emission comprises ~ 8% of the total 4- 
20 keV flux. Interpreting this emission as due to Compton scattering of 
relativistic electrons (which produce the known extended radio emission) 
by the cosmic microwave background radiation, we determine that the 
mean, volume- averaged magnetic field in the central region of Coma is 
5 -0.1- 0.3 /iG. 



Subject headings: Calaxies: clusters: general — galaxies: clusters: in- 
dividual (Coma) — galaxies: magnetic fields — radiation mechanisms: 
non-thermal 
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1. Introduction 

X-ray spectra of clusters of galaxies have 
long been expected to show structure be- 
yond that of a single temperature thermal 
model, mainly due to non-isothermality of in- 
tracluster (IC) gas in the outer cluster region. 
In addition, non-thermal (NT) X-ray emis- 
sion in clusters was predicted {e.g., , Rephaeli 
1977) from Compton scattering of relativis- 
tic electrons by the Cosmic Microwave Back- 
ground (CMB) radiation. There is at least 
some observational evidence for radial varia- 
tion of the gas temperature in a few clusters 
{e.g.,, Markevitch 1996, Honda et a/. 1996, 
Donnelly et a/. 1999, Watanabe et a/. 1999). 
In the Coma cluster, recent XMM measure- 
ments indicate (Arnaud et al. 2001) that the 
temperature is remarkably constant within 
the central region where temperature varia- 
tion was previously deduced from ASCA mea- 
surements. After a long search (for a re- 
cent review, see Rephaeli 2001), NT X-ray 
emission seems to have finally been measured 
in Coma (Rephaeli, Gruber & Blanco 1999, 
hereafter RGB, Fusco-Femiano et a/. 1999), 
A2256 (Fusco-Femiano 2000), A2319 (Gruber 
& Rephaeli 2002), and perhaps also in A2199 
(Kaastra et a/. 2000). 

Appreciable deviation from isothermality 
may have significant impact on modeling the 
structure and evolution of IC gas, and on use 
of the gas as a probe to determine the to- 
tal cluster mass (assuming hydrostatic equi- 
librium). The exact gas density and tem- 
perature profiles are also very much needed 
in analysis of measurements of the Sunyaev- 
Zeldovich (S-Z) effect and its use as a cosmo- 
logical probe. There clearly is strong moti- 
vation for a more realistic characterization of 
cluster X-ray spectra for an improved descrip- 
tion of IC gas, and the study of NT phenom- 



ena in clusters. 

We have previously analyzed ~ 87 and 
~ 160 ks RXTE measurements of the Coma 
cluster and A2319, respectively, in order to 
search for NT emission from these clusters 
which have well documented extended regions 
of radio emission. Analyses of these measure- 
ments yielded strong evidence for a second 
spectral component in both clusters. While 
the second component could possibly indicate 
a temperature variation across the cluster, we 
have concluded (RGB, Gruber & Rephaeli 
2002) that the deduced spectral parameters 
are more consistent with power-law emission. 
NT emission in Coma seems to have been de- 
tected directly - in the 25-80 keV range - by 
BeppoSAX (Fusco-Femiano et a/. 1999). This 
provided further impetus to propose a longer 
observation of this cluster with RXTE. Here 
we briefly report the results from a joint anal- 
ysis of these and the previous RXTE measure- 
ments, with a total integration time of ~ 264 
ks, and discuss some of their direct implica- 
tions. 

2. Observations and Data Reduction 

Coma was observed with the Proportional 
Counter Array (PCA) and the High Energy 
X-ray Timing Experiment (HEXTE) on RXTE 
during 58 separate pointings totaling approx- 
imately 87 ks in June 10 - 22, and July 15, 
1996, and during 50 additional pointings to- 
taling nearly 177 ks in November 24 - Decem- 
ber 15, 2000. Spectral results from the earher 
observation were reported by RGB. 

PCA data were collected in the 'Standard 
2' spectral mode, which consists of a 129- 
channel count spectrum nominally spanning 
2 to 1000 keV with 16 second time resolu- 
tion. Data from the two independent HEXTE 
clusters of detectors were taken in event-by- 
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event mode, which were subsequently accu- 
mulated into 256-channel spectra spanning 
17-250 keV. To subtract the background, each 
HEXTE cluster was commanded to beam- 
switch every 16 s or 32 s between on-source 
and alternate off-source positions 1.5° on ei- 
ther side. 

Standard screening criteria were applied 
to the data segments (Earth elevation angle, 
spacecraft pointing, avoidance of the South 
Atlantic Anomaly, times of geomagnetic ac- 
tivity), resulting in a net exposure time for 
PCA of 87328 s during the 1996 observations 
and 176864 s during the second observation 
set. During the first set of observations PCA 
detectors 0, 1 and 2 were always on, but de- 
tectors 3 and 4 were enabled for less than 
half of the observation set, and were ignored 
in the interest of a more uniform data set. 
The (year) 2000 observations were made with 
PCA detectors and 2 only, but the detec- 
tor spectra could not be used because of 
the loss of the guard propane layer, result- 
ing in higher and uncertain background. The 
detector 2 gain was about 10% higher in the 
later observations. The HEXTE data, by con- 
trast, formed a perfectly uniform data set, 
with no changes of number of detectors or 
gain. HEXTE net times were shorter than 
PCA by nearly a factor of four: half the time 
was used for background measurement, and 
nearly another factor of two was lost due to 
electronic dead time caused by cosmic rays. 

The PCA background was estimated with 
the 'L7/240' faint source model provided by 
the instrument team. No significant emission 
from Coma was detected above ~ 40 — 50 
keV. The spectral data above 40 keV were 
used to determine a correction of ~0.5% to 
the background estimate. No such tuning was 
required for the HEXTE background, which 



is measured, and has been determined (Mac- 
Donald 2000) to be accurate to within a few 
hundredths of a count /s in long exposures. 

3. Spectral Analysis 

At about 40 count/s per PCA detector, the 
observed counting rate from Coma was well 
above the background rate of 13 count/s per 
detector, thus the cluster was easily detected 
and possible errors of background estimation 
are unlikely to compromise spectral analy- 
sis. The HEXTE net rate of somewhat less 
than one count/s per cluster was well below 
the background rate of 80 counts/s per clus- 
ter but still large compared to the error limit 
of roughly 0.03 counts/s determined by Mac- 
Donald (2000). Inspection of the screened 
light curves for PCA and HEXTE revealed no 
significant variation, as expected for a cluster 
of galaxies. Accordingly, we co-added all the 
selected PCA and HEXTE data to form net 
spectra for analysis. Spectra from the indi- 
vidual PCA detectors were also coadded, and 
the small differences in gain were accounted 
for in the generation of the energy response 
matrices, following procedures prescribed by 
the PCA analysis team. 

A small energy dependence for systematic 
errors (e.g., Gruber et al 2001, Wilms et al. 
1999) averaging about 0.3%, was applied to 
the PCA spectrum. Additionally, PCA spec- 
tral channels below 4keV and above 22keV 
were excluded because of sensitivity to arti- 
facts in the background model, as well as the 
rapidly dechning effective area of the PCA 
outside these bounds. The HEXTE data 
were restricted to the energy range 19-80 keV 
for similar reasons, resulting in source and 
background counting rates of 0.75 and 84.1 
count/s respectively (cluster A), and 0.49 and 
58.1 count/s (cluster B). 
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Repeating the procedure in RGB, we fit the 
joint PGA and HEXTE spectra summed over 
all observations to three simple spectral mod- 
els: a Raymond-Smith (R-S) thermal plasma 
emission model, R-S plus a power law model, 
and two R-S models at different tempera- 
tures. In all three cases most of the observed 
flux is in a primary ~ 8 keV R-S component. 
Best-fit parameters and 90% confidence inter- 
vals are hsted in Table 1. The best-fit tem- 
perature in the single isothermal gas model is 
7.90 ± 0.03 keV, in good agreement with the 
range determined from ASGA measurements 
(Honda et al. 1996), and only ~ 4% less than 
the value of 8.2 keV deduced recently from 
XMM measurements of the central it/ 'region 
of Goma (Arnaud et a/. 2001). This region 
is much smaller than the RXTE ~ 1° field 
of view. The observed ~ 6.7 keV Fe XXV 
Kq, line yields an abundance of 0.195 ± 0.008 
(in solar units), quite consistent with previ- 
ously determined values (> 0.2). No cold ab- 
sorption was measurable, and given the 4 keV 
PGA threshold, none was expected. 

The poor quality of the fit to a single 
isothermal model is apparent from the fact 
that = 76.8 for 45 degrees of freedom. 
Residuals have a high-low-high pattern which 
signals the need for another smooth spectral 
component. When a second thermal compo- 
nent is added, best-fit parameters are kTi ~ 
7.5 keV, and a very high kT2 ~ 37.1 keV, 
with the second component accounting for a 
modest fraction, ~ 6% of the total fiux. For 
this fit, = 48.1, lower by 28.8 than the 
value obtained in the single R-S model. The 
F-test probabihty of the second component 
is 0.9992 for two additional degrees of free- 
dom. The range of solutions is very large, 
however, because the problem is numerically 
highly degenerate. If we consider the possible 



sets of (kTi, kT-i), the contour + 4.6 de- 
fines the joint 90% probability contour. The 
lowest temperature combination permitted in 
this range is kTi = 5.48 keV and kT2 = 9.0 
keV, with fractions of the 4-20 keV flux, re- 
spectively, of 24% and 76%. The highest tem- 
perature combination has kTi ~ 7.5 keV and 
kT2 unbounded. Equal 4-20 keV flux contri- 
butions are obtained with the values 6.5 and 
10.5 keV for the two temperatures. 

A somewhat better flt than the two-temperature 
model, Ax^ = 30.7 (with respect to the sin- 
gle R-S model), is obtained when the second 
component is a power-law, with a photon in- 
dex of 2.1 ± 0.5 (90% confldence). The 4-20 
keV flux of the power-law component is 8% of 
the total. This component comprises most of 
the flux only at energies > 60 keV. The iron 
abundance changes negligibly from the value 
obtained with the isothermal flt. The count 
rate from the combined thermal and power- 
law emissions, and that of just the power-law 
component, are shown in Figure 1, together 
with the measurements. 

Only very limited comparisons can be made 
with other satellite measurements due to dif- 
ferences in spectral bands and fleld of views. 
Using archival ASGA GIS spectra of the cen- 
tral 20' diameter region of Coma, we flnd 
that these are well flt by a single tempera- 
ture Raymond-Smith model, with x^ = 664 
for 713 degrees of freedom. However, fltting 
an extra power law component with number 
index 2.3 reduces x'^ by 12, and thus is sig- 
niflcant at about the 3a level. The 3-10 keV 
band common to both ASGA and RXTE is 
most useful for comparison: the ASGA power 
law represents 5 ± 2 percent of the total emis- 
sion, consistent at 68% with the correspond- 
ing RXTE value of 7%. 

EUV emission from the Goma cluster was 
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Table 1: Results of the spectral analysis 



Parameter 


single R-S two R-S R-S+power-law 


xVdof 
kTi (keV) 

Normalization" 
Abundance^ 
m (keV) 

Normalization" 
/,(5 keV) (cm-2s-i) 
Photon index 


76.8/45 48.1/43 46.1/43 
7.90 ±0.03 7.47 ±0.22 7.67±0.1 
0.348±0.002 0.345±0.020 0.334±0.016 
0.195±0.008 0.192±0.010 0.202±0.017 
37 ±29 
0.011±0.011 

(3.3±1.0) X 10-^ 
2.07±0.46 



Notes: 

All quoted errors are at the 90% confidence level, 
"e.m. = Raymond-Smith emission measure in units of 10"^'^ J NgNndV j AnD^, where D is the luminosity distance 

and Ne, Nh are the total number of electrons and protons, respectively. 
''Abundance is expressed relative to solar values. 



measured with the EUVE by Bowyer & Berghofer 
(1998). The emission is said to be at a level 
significantly higher than what is expected 
from the main X-ray emitting gas, and pre- 
sumably cannot be explained simply by in- 
voking a second, colder, T ~ 10^ K, gas 
component. A similar conclusion was re- 
cently reached by Durret et al. (2002) from re- 
analysis of EUV and ROSAT measurements. 

A more direct comparison of our spectral 
results can be made with BeppoSAX obser- 
vations of Coma with the MEGS and PDS 
experiments. Analysis of the PDS measure- 
ments led Fusco-Femiano et al. (1999) to con- 
clude that a significant (at the 4. 60" level) NT 
flux, 2.2x 10~^^ ergs cm"^ s~^, was detected at 
energies 20—80 keV. While this range includes 
the values we inferred from the RXTE mea- 
surements, it should be noted that an accu- 
rate spectral analysis of the combined MEGS 
and PDS data is not possible because of the 
very different fields of view of these experi- 
ments. 



4. Discussion 

The consistent results from the analysis of 
the second (year 2000) and combined (years 
1996 & 2000) datasets further substantiate 
the reality of the detection of a second compo- 
nent in the Goma spectrum. If thermal, this 
component could help determine the thermal 
structure of IG gas, which in turn would have 
important implications for the use of clus- 
ters as cosmological probes in general, and 
the most extensively researched nearby rich 
cluster in particular. If this component is 
NT emission from a population of relativis- 
tic electrons, then a new dimension for the 
study of NT IG phenomena will have been 
opened. Since the combined ~ 264 ks RXTE 
observations have not yielded direct detection 
of power-law emission at energies > 30 keV, 
we have to invoke other observational and 
theoretical considerations in order to identify 
more uniquely the origin of the extra emission 
we have deduced. 

It is unrealistic to expect that IG gas is 
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Fig. 1.— The RXTE (photon) spectrum of 
the Coma cluster and folded Raymond-Smith 
(kT ~ 7.67), and power-law (index = 2.3) 
models. HEXTE data points are marked with 
circles and 68% error bars. The total fitted 
spectrum is shown with a histogram, while 
the lower histogram shows the power-law por- 
tion of the best fit. The quality of the fit is 
demonstrated in the lower panel, which dis- 
plays the observed difference normalized to 
the standard error of the data point. 



fully isothermal outside the central region of 
the cluster. A more likely behavior is at least 
some decrease of the temperature outside a 
central ~ 2 — 3 core radii region. Indeed, pre- 
vious ASCA measurements of Coma (Honda 
et al. 1996) and other clusters (e.g., Marke- 
vitch 1996) seem to have shown some devia- 
tions from isothermality. In Coma, the Honda 
et al. (1996) analysis of measurements from 14 
different pointings of an area extending ~ 1° 
from the center indicated that the IC gas tem- 
perature varied by ±50% - with respect to 
the overall mean value of ~ 8 keV - in two 
azimuthal regions 40' from the center. How- 
ever, these two regions cover only a small part 
of the projected area of the cluster, and since 
the temperature is higher than the mean in 
one region, while it is lower in the other, the 
overall change of the spectral flux, as com- 
pared with that from an isothermal gas at the 
mean cluster temperature, is negligible. More 
recent results from high spectral and spatial 
resolution measurements of Coma with XMM 
it was concluded that the temperature is con- 
stant in the central ~ 10' radial region, with 
a best fit value of 8.2 ±0.1 keV (Arnaud et 
al. 2001). While the fits to the XMM data do 
not seem to have included two-temperature 
models, it is clear that an appreciable emis- 
sion at a significantly different temperature 
than this mean value would have been de- 
duced from the XMM measurements through 
a larger variance, if not in the form of a sys- 
tematic temperature gradient. 

The results of our RXTE analysis yield a 
statistically most probable temperature com- 
bination with a prohibitively high value for 
the second component, kT2 ~ 37.1. At its 
lowest boundary, the 90% contour region in 
the {kTi, kT2) plane does include the more ac- 
ceptable values kTi = 5.5 keV, and kT2 = 9, 
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but this only if the respective 4-20 keV frac- 
tional fluxes of these two components are 24% 
and 76%. It is quite unlikely that about a 
quarter of the flux could come from a compo- 
nent with a significantly lower temperature 
than the mean value deduced by virtually 
all previous X-ray satellites. In particular, 
such a component would have been detected 
in the high spatially resolved measurements 
with ROSAT and XMM. 

To assess the possibihty that a two-temperature 
gas model is just a simplified representation 
of a more realistic continuous temperature 
distribution, we have repeated the following 
simple procedure we employed in our anal- 
ysis of the first RXTE observations (RGB): 
Assuming a polytropic gas temperature pro- 
file of the form T(r) oc n(r)'''~"^, with the 
familiar /? density profile for the gas den- 
sity, n(r) oc (1 -l- r^/r^)~^^/^, where Tc is 
the core radius, we calculated the integrated 
flux and the mean emissivity-weighted tem- 
peratures as functions of 7, f3, and r. These 
quantities were then calculated in the regions 
[0, r] and [r, Rq] by convolving over the tri- 
angular response of the PCA with Rq ~ 58'. 
Prom ROSAT observations, Tc ~ 10.0', and 
p ~ 0.70d=0.05 (Mohr et al. 1999). We sought 
the range of values of r, /5, and 7 for which the 
two mean emissivity-weighted temperatures 
and respective fluxes from these regions are 
closest to the values deduced from our spec- 
tral analysis in Section 3. The results of these 
calculations indicate that for 0.5 < (3 < 0.9 
and 1 < 7 < 5/3, there is no acceptable poly- 
tropic configuration that matches the obser- 
vationally deduced values of the temperatures 
and fractional fiuxes. For low values of 7 the 
temperature gradient is too shallow, while for 
high values the implied central temperature 
is unrealistically high. This simple plausibil- 



ity check suggests that the two thermal com- 
ponents model is somewhat inconsistent with 
the RXTE results. However, the gas distribu- 
tion may be more complicated than consid- 
ered here, so that a temperature structure as 
implied here cannot be altogether ruled out. 

Of particular interest is the somewhat more 
likely possibility that the secondary spectral 
component is NT. Since emission from an 
AGN in the FOV is not likely (see details 
in RGB), it is natural to consider that this 
emission is due to Compton scattering of rel- 
ativistic electrons whose presence in Coma 
is directly inferred from many measurements 
of spatially extended region of radio emission 
{e.g., Kim et al. 1990, Giovannini et al. 1993). 
Prom the measured radio spectral index, 1.34± 
0.1, it readily follows that the predicted power- 
law (photon) fiux from Compton scattering 
of these electrons by the CMB has an in- 
dex 2.34 ± 0.1, a value which is quite con- 
sistent with what we have inferred, 2.1 ± 0.5 
(all errors arc at 90% confidence). With 
the measured mean radio fiux of 0.72 ± 0.21 
Jy at 1 GHz, the power-law X-ray fiux de- 
duced here, and the assumption that the spa- 
tial factors in the theoretical expressions for 
the two fiuxes are roughly equal, we can eas- 
ily compute (see more details in RGB) the 
mean volume-averaged value of the magnetic 
field, Brx- Taking into account the full 90% 
range of values of the radio flux, radio in- 
dex, and the power-law X-ray fiux, we get 
B.f.rj. ~ 0.1 — 0.3/iG. This range of values for 
B.f;j. is consistent with our previous estimate 
(RGB), and the range (0.14 - 0.25 /xG) de- 
duced by Pusco-Pemiano et al. (1999). 

Since we have assumed that the spatial fac- 
tors in the theoretical expressions for the ra- 
dio and NT X-ray fiuxes are roughly equal, 
it follows that the mean value of the deduced 
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magnetic field is independent of the source 

size and distance. To determine the rela- 
tivistic energy density we do have to specify 
the radius of the emitting region. Scahng to 
the observed radius of the diffuse radio emis- 
sion, R ~ 20', and integrating the electron 
energy distribution over energies in the ob- 
served radio and X-ray bands, we obtain pe — 
(8 ±3) X W-^\R/20')-^ ergcm-3; a distance 
of 139 Mpc (with Ho = 50 kms"^ Mpc"^) was 
used. Based on the high Galactic proton to 
electron energy density ratio of cosmic rays, 
it can be conjectured that the energetic pro- 
ton energy density is considerably higher than 
this value. 

The strength of IC magnetic field can also 
be estimated from Faraday rotation measure- 
ments of background radio sources seen through 
clusters, yielding a different mean field value, 
Bfr- Analyses of such measurements usu- 
ally yield field values that are a few /iG 
(see, e.g.,, Clarke, Kronberg, and Bohringer 
2001, and the review by Carilli & Taylor 
2002). Clearly, the mean strength of IC 
fields has direct implications on the range of 
electron energies that are deduced from ra- 
dio measurements, and therefore on the elec- 
tron (synchrotron and Compton) energy loss 
times. Higher electron energies imply shorter 
energy loss times, with possibly important 
ramifications for relativistic electron models 
{e.g., , Rephaeli 1979, Sarazin 1999, Ensslin et 
al. 1999, Brunetti et al. 2001, Petrosian 2001). 

Much has been written about the apparent 
discrepancy between deduced values of -B^x 
and Bfr- Indeed, it is sometimes claimed that 
this discrepancy makes Compton interpreta- 
tion of cluster power-law X-ray emission un- 
tenable. However, B^x and Bfr are actually 
quite different measures of the field: Whereas 
the former is essentially a volume average of 



the relativistic electron density and (roughly) 
the square of the field, the latter is an average 
of the product of the line of sight component 
of the field and gas density. All these quan- 
tities vary considerably across the cluster; in 
addition, the field is very hkely tangled, with 
a wide range of coherence scales which can 
only be roughly estimated. These make the 
determination of the field by both methods 
considerably uncertain. Thus, the unsatisfac- 
tory observational status (stemming mainly 
from lack of spatial information) and the in- 
trinsic difi'erence between B^x and Bf^, make 
it clear that these two measures of the field 
cannot be simply compared. Even ignoring 
the large observational and systematic uncer- 
tainties, the different spatial dependences of 
the fields, relativistic electron density, and 
thermal electron density, already imply that 
Brx and Bfr will in general be quite different. 
This was specifically shown by Goldshmidt & 
Rephaeli (1993) in the context of reasonable 
assumptions for the field morphology, and the 
known range of IC gas density profiles. It 
was found that Brx is typically expected to 
be smaller than Bfr. Various statistical and 
physical uncertainties in the Faraday rotation 
measurements, and their impact on deduced 
values of IC fields, were investigated recently 
by Newman, Newman & Rephaeh (2002); 
their findings further strengthen the assess- 
ment that a simple minded comparison of val- 
ues of Brx and Bfr is meaningless, and that 
it is quite premature to draw definite conclu- 
sions from the apparent discrepancy between 
values deduced by these very different meth- 
ods to measure IC magnetic fields. 

As we have mentioned in the previous sec- 
tion, the excess EUV emission in Coma could 
also be NT, and based on the similar mor- 
phologies of the EUV emission and low fre- 
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quency radio emission, Bowyer & Berghofer 
(1998) interpreted this emission as Compton 
scattering of the CMB by a population of low 
energy electrons. They adopted a value for 
the power-law index which is somewhat lower 
than the value used here, but deduced a sim- 
ilar value (~ 0.2 /iG) for the mean magnetic 
field. More recently, Tsay, Hwang & Bowyer 
(2002) have explored whether a Compton ori- 
gin for the observed EUV excess can be main- 
tained even if the field is as high (few /^G) 
as is currently deduced from Faraday rota- 
tion measurements. They conclude that this 
is possible within a limited class of lower en- 
ergy (~ 100 MeV) NT electron models, and 
that in this different (second) popu- 

lation of relativistic electrons is required to 
explain the measurements of NT X-ray emis- 
sion by RXTE and BeppoSAX. Note that sig- 
nificant IC density of sub-relativistic electrons 
could in principle also produce high energy X- 
ray emission by NT bremsstrahlung (Kaastra 
et al. 1998, Sarazin & Kempner 2000). How- 
ever, the properly normalized contribution of 
such electrons to the power-law emission de- 
duced here from the RXTE measurements is 
too small (Shimon & Rephaeli 2002) to affect 
our estimated value of the magnetic field. 

Further evidence for the NT nature of the 
second spectral component in Coma could 
possibly come from the scheduled 500 ks ob- 
servation of this cluster with with IBIS imager 
aboard the INTEGRAL satellite. The mod- 
erate ~ 12' spatial resolution of IBIS can po- 
tentially yield crucial information about the 
location and size of high energy NT X-ray 



emission. 



We are grateful to the referee for useful 
comments made on an earlier version of the 
paper. 
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